Signatures of evanescent mode transport in graphene 
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We calculate the shot noise generated by evanescent modes in graphene for several experimental 
setups. For two impurity-free graphene strips kept at the Dirac point by gate potentials, separated 
by a long highly doped region, we find that the Fano factor takes the universal value F = 1/4. For 
a large superlattice consisting of many strips gated to the Dirac point interspersed among doped 
regions, we find F = l/(81n2). These results differ from the value F = 1/3 predicted for a disordered 
metal, providing an unambiguous experimental signature of evanescent mode transport in graphene. 

PACS numbers: 72.80.Vp, 73.23.Ad, 73.50.Td, 73.63.-b 



I. INTRODUCTION 

The nonequilibrium current fluctuations, or shot noise, 
in ballistic graphene^ - — have received much attention 
since the seminal paper by Tworzydlo et al£> In that 
work, it was shown that shot noise can be generated even 
in a completely impurity-free sheet of graphene. This 
result is surprising at first sight, as the shot noise van- 
ishes in conductors without electron scattering^— The 
unanticipated noise is caused by evanescent (exponen- 
tially damped) waves that backscatter electrons, even in 
clean graphene. 

It has been shown by Tworzydlo et aL— that in a clean 
sheet of graphene at its Dirac point, that is, at zero chem- 
ical potential, the shot noise is not only nonzero, but, 
moreover, it has universal characteristics. The shot noise 
normalized by the mean current and expressed in units of 
the electron charge e - this quantity is commonly referred 
to as the Fano factor F - takes the value F = 1/3. This 
prediction has generated much theoretical interest in the 



topic, 
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and considerable experimental activity. 
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So far, the prediction of Tworzydlo et ai£- has been 
tested by two experiments, as reported in Refs.[l8landlT9l. 
and a Fano factor close to F = 1/3 was measured. How- 
ever, the interpretation of these experiments is ambigu- 
ous, as the universal shot noise value for evanescent mode 
transport at the Dirac point of graphene^ is identical to 
the one expected for conventional disordered conductors 
in the diffusive regime^ It has been confirmed in nu- 
merous regimes, both witbJ^ and without - — electron- 
electron interactions, that a Fano factor close to F = 1/3 
is also expected for disordered graphene, independent of 
the chemical potential. Therefore, the measured Fano 
factor F = 1/3 could be due to either evanescent mode 
transport or impurities in the measured samples. 

In the experiment of Danneau et al.,^ an additional 
signature of evanescent mode transport was observed; 
namely a strong dependence of the measured Fano factor 
on the chemical potential. Such dependence on the chem- 
ical potential is expected for clean graphenei 5 ' 15 ! 16 As the 
chemical potential departs from the Dirac point, more 
and more evanescent waves become propagating and 
cease to backscatter electrons, decreasing the Fano factor. 
The observed dependence of F on the chemical potential 



does not occur in generic diffusive conductors ji^r— mak- 
ing it a more distinctive signature of evanescent mode 
transport. Nevertheless, other scenarios are also con- 
sistent with the doping dependence of F reported by 
Danneau et al.^ such as energy- dependent scattering. 
Additional experimental signatures of transport through 
evanescent modes in graphene are therefore desirable. 

In this article, we propose experimental geometries for 
which the transport through evanescent modes at the 
Dirac point of graphene has unambiguous signatures. We 
first study a sheet of graphene subject to gate potentials 
that induce two strips of graphene with chemical poten- 
tial at the Dirac point, separated by a highly doped re- 
gion. We show that the Fano factor takes a universal 
value in this geometry, as in that of Ref. S F = 1/4. 
This value is different from that of diffusive conductors, 
and it provides an unambiguous signature of evanescent 
modes. Graphene devices with multiple potential steps, 
as required for the proposed test of evanescent transport, 
have been implemented experimentally^ 

Similar results are obtained for longer cascades of p > 2 
strips of graphene with evanescent transport. In par- 
ticular, we take the limit p — > oo of a long graphene 
superlattice with a piecewise constant potential. Such 
superlattices have recently received much attention as a 
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way to engineer the bandstructure of graphem 
the point of creating new Dirac cones in the electronic 
spectrum^ - — Under certain conditions, we find for such 
superlattices yet another universal value of the Fano fac- 
tor: F = l/(81n2). 



This paper is organized as follows: In section |TT1 we 
analytically derive the universal Fano factor F = 1/4 for 
two evanescent regions in series. We study numerically 
the departures of F from this value as asymmetries are 
introduced into the setup, an additional experimental sig- 
nature of evanescent mode transport. In section [TTTl we 
first generalize our approach to cascades of an arbitrary 
number p of evanescent regions in series before taking 
the long superlattice limit p — > oo. We conclude with a 
summary in section IIVI 
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FIG. 1. Schematic of a graphene strip of width W consisting 
of a highly doped region (gray rectangle) of length A at voltage 
Vd separated by weakly doped regions of length L and R. 
Separate gate leads (not shown) control the voltages Vl and 
Vr in the weakly doped regions. The system is contacted at 
either end by ideal leads (not shown) at the lead potential Vg. 
The electrostatic potential in the strip is plotted above the 
schematic as a function of x. 



II. EVANESCENT MODE TRANSPORT 
THROUGH TWO EVANESCENT STRIPS IN 
SERIES 

We analyze the transport through a graphene sheet of 
width W subject to gates that allow one to tune two 
strips of lengths L and R to a chemical potential close to 
the Dirac point. We denote the electrostatic potentials in 
those strips by Vl and Vr, respectively. The two strips 
are separated by a highly doped region of length A at 
electric potential Vd (see Fig. [T]) . 

We model the leads that contact the sample by highly 
doped regions of graphene to the left and to the right 
of the sample at voltage Ve, as in the calculations of 
Tworzydlo et al£ It has been confirmed in density func- 
tional calculations that this model correctly describes 
transport into certain types of contacts^ In units where 
h = 1 , our model Hamiltonian takes the form 

Hj — VCTj ■ p + V(x), (1) 

where <r 7 = (^ya x ,a v ) is a vector of Pauli matrices with 
the valley index 7 = ±1, p is the electron momentum and 
v is the electron velocity. The potential V(x) takes the 
value Vd in the highly doped regions, Vg in the leads and 
the gate voltages Vl and Vr in the regions that can be 
tuned to the Dirac point (see Fig.[T]). We assume the mi- 
croscopic potential, which is represented by the potential 
V(x) in our long-wavelength theory, to be smooth on the 
lattice scale, so that it does not scatter between valleys. 

The plane wave solutions of the Dirac equation have 
energy e = Vd ± v \Jk 2 d + q 2 in the highly doped regions 



and e — Vl(r) ± v yk L (R) + In m tne ( r i§ nt ) gated 
regions, where the ± sign refers to the conduction and va- 
lence bands, respectively. We choose V^(<q 3> v/L,v/R 
such that, at the Fermi level, all relevant modes are prop- 
agating (real k^)) in the leads and the highly doped 
regions. In the weakly doped regions, both propagating 
(real kL(R)) and evanescent (imaginary kL(R)) modes can 
occur. 

The transverse wavenumbers q n depend on the bound- 
ary conditions at y = and y = W . We consider a class 
of boundary conditions that do not couple longitudinal 
and transverse wavenumbers^ In this case, the transverse 
wavenumber q n is a good quantum number, and modes 
are not mixed at the weakly-to-highly doped region in- 
terfaces. For the infinite mass, metallic armchair edge, or 
semiconducting armchair edge boundary conditions, the 
transverse wavenumber is given by q n = (n + a)ir /W , 
where^ 

{1/2 infinite mass 
metallic armchair 
1/3 semiconducting armchair. 

In the continuum limit W 3> L,R, which we take for 
the remainder of this article, one expects the transport 
properties to be independent of the boundary conditions. 
Accordingly, a drops out of the calculation. 



A. The universal limit 

We first consider the situation in which both weakly 
doped regions of the sample are tuned to the Dirac point, 
Vl = Vr — 0. Matching modes at the interfaces between 
segments of differing potential results in the transmission 
amplitude 

t n = [e lkdX sinh(n£) sinh(n^) + 

+e -ik d X cosn (n£) cosh(n-ft)] _1 (2) 

for mode n at the Fermi level, which we choose to be at 
energy e-p = 0. The dimensionless "lengths" C and 1Z are 
given by £ = nL/W and K = ttR/W. The Fano factor 
is found from the transmission probabilities T n = \t n \ 2 
asi^ 



Sri -^n 

In the continuum limit, W 3> L, R 7 the sums become 
integrals over the mode index, and the Fano factor is 
given by F = 1 — I2 / 1\ with 

J i = Jo°° T ndn, 

I 2 =J~T*dn. (4) 

In Fig. [21 the transmission probability determined from 
Eq. @ is integrated numerically in order to obtain the 
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FIG. 2. Fano factor F for two graphene strips as in Fig.[T]for 
L = R and chemical potential at the Dirac point, as a function 
of the length of the highly doped region in units of \nd\ L 2 . 
The curve is calculated from Eq. Q in the limit N — > oo. The 
Fano factor approaches 0.33 .. . as A — > 0, in agreement with 
the calculations of Tworzydlo et al^ For large values of the 
length of the highly doped region, the Fano factor approaches 
0.25 . . ., in agreement with the analytic calculations presented 
in the text. 



Fano factor for a symmetric system (L = R) as a func- 
tion of the thickness A of the central, highly doped re- 
gion. As A — >• 0, the Fano factor approaches 1/3. This 
result agrees with the calculations of Tworzydlo et al., 5 
as this limit corresponds to transmission through a single 
graphene strip at the Dirac point. In contrast, the Fano 
factor approaches 0.25 ... in the limit A ;§> \nd\ L 2 , where 
Kd = V d /v. 

This numerical result suggests that the Fano factor 
may be accessible analytically in the limit A 3> \nd \ L 2 . 
Indeed, one finds that the integrals in Eq. ([J| can be done 
analytically for A 3> \n d \ L 2 . Consider first I\. We write 
the transmission probability determined from Eq. ([2]) in 
the form 



T„ 



1 



1 



a(n) l + j3(ri) cos[2k d (n)X] , 



(5) 



where 



a(n) — cosh 2 (n£) cosh 2 (n7?.) + sinh 2 (n£) sinh 2 (n7?.), 
/3(n) = 2cosh(n£) cosh(n7?.) sinh(n£) s'mh(nlZ) / a(n) 

k d (n) = ^K 2 d - n 2 n 2 /W 2 . 



(6) 



The transmission probability decays exponentially with 
n(£ + 1Z), and the integrals in Eq. (j3J are thus cut-off at 
n ~ max{l/£, 1/72.}. 

The key observation is that for A ^> \Kd \ L 2 , the cosine 
function in Eq. ([5]) oscillates rapidly on the scale of the 
exponential decay of T n (recall that we assume \n d \ L 3> 
1). In the limit A/ \Kd\ L 2 — > oo, therefore, a(n) and /3(n) 



An(n ) 



KdW 2 
Xttuq 



(7) 



of the cos [2kd(n)X] function around a given index n . 
Consequently, I\, restricted to an interval of length 
An(no) around no, becomes the integral of the func- 
tion [1 + (3(n) cos (5 — unjf 1 , which can be done analyt- 
ically. Here, S and uj are found by linearizing 2kd{n)\ = 
6{n ) - Lu(n )n with 6 (no) = 2n d \(l + n 2 a T: 2 /2n 2 d W 2 ) 
and u(no) = 2tt 2 \n / HdW 2 . In the limit \/\nd\L 2 -> 
oo, the periods An(n ) are short compared to the de- 
cay scale of the transmission probability set by L + 1Z, 
and the requisite sum over all periods becomes a second 
integral. Noting that \(3(n )\ < 1, this results in 



h 



dnn 



(8) 



where 



7±( n o) = cosh 2 (n £) cosh 2 (n 7?.) ± 
± sinh 2 (n £) sinh 2 (n 7?.). 



(9) 



The integral I2 of Eq. can be done in the same man- 
ner, with the result 



h 



7+Oo) 
7-(»o) 



driQ. 



(10) 



In the symmetric case, L = R, the integrals in Eq. 
([8]) and (JTO]) can also be evaluated analytically, result- 
ing in the same Fano factor as for a ballistic quantum 
dotj2i F = 1/4, in accordance with the numerical results 
shown in Fig. [5] Evanescent transport in this geometry 
thus has an unambiguous signature, with a Fano factor 
that differs from the one in a disordered sample. For 
the asymmetric case, we calculate the Fano factor nu- 
merically as a function of R/L, with the results plotted 
in Fig. [3l The Fano factor approaches 1/3 in the limit 
R/L — > 00, which again corresponds to a single graphene 
region at the Dirac point, as considered by Tworzydlo et 

The dependence of the Fano factor on the ratio R/L 
shown in Fig. [3] is a more distinctive signature of evanes- 
cent wave transport than the value F— l/Aa,tL = R 
alone. In a typical experiment, however, the lengths of 
the gated regions cannot be changed easily. Alterna- 
tively, the gate voltages Vl and Vr can be controlled. 
We discuss the dependence of F on those gate voltages 
in the following subsection. 



B. Voltage induced signatures of evanescent mode 
transport 

Here, we obtain the dependence of the Fano factor on 
the gate voltage Vr in the symmetric configuration, L = 
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FIG. 3. Length ratio dependence of the Fano factor for 
two graphene strips as in Fig. [T] with chemical potential at 
the Dirac point. The curve is calculated from the integrals 
given in Eq. Q and Eq. (| 10p in the continuum limit with 
A ^> \Kd\ L 2 ■ The Fano factor shows a minimum of 0.25 . . . 
for the symmetric case {R/L = 1) and tends toward 0.33 . . . 
as R/L -> or R/L — > oo, in agreement with the results of 
Tworzydlo et al^ 



R, (cf. Fig. [IJ with the left gated region at the Dirac 
point, Vl=0. The transmission amplitude as a function 
of Ur takes the form 



] tR [g ik d \ CQS ^ n £ CQS ^ R l_ 



-ikr sin UrL) + e d — smb. nC sin kjiL 



(11) 



where n 2 R(d) 



k R(d) +1n with K Bid) = V R{d) /v. We 



sider the limit A/ \nd\ L 2 —> oo of the resulting expression 
for the Fano factor, following the methods of the previ- 
ous section. In this case, we perform the final integration, 
corresponding to the sum over periods An(no), numeri- 
cally, with the result plotted in Fig. |4j The Fano factor 
is 1/4 when the gate voltage Vr is at the Dirac point, as 
previously calculated, and it approaches 1 /3 in the limit 
Vr 3> v/L. The crossover between F = 1/3 and F = 1/4 
within a gate voltage interval of AV"l ~ v/L serves as an- 
other distinctive signature of evanescent mode transport. 



III. EVANESCENT MODE TRANSPORT IN A 
GRAPHENE SUPERLATTICE 



0.34r 



0.32 



0.28 



0.26 



0.24 L 



4 6 

V R L/v 



10 



FIG. 4. The Fano factor for two graphene strips as in Fig. [T] 
for L = R and A 3> \kJ\ L 2 as a function of gate voltage Vr 
in the right graphene region, measured in units of v/L. Gate 
leads fix the chemical potential in the left evanescent region 
to the Dirac point (Vl ~ 0). The Fano factor is 0.25 . . . 
when Vr is at the Dirac point, as previously calculated, and 
it approaches 0.33 . . . as Vr 2> v/L. 
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FIG. 5. Schematic of a segment of a graphene superlattice of 
width W , consisting of evanescent graphene regions of length 
L separated by doped regions (gray rectangles) of length A 
at the potential Vd- Separate gate leads (not shown) fix the 
voltages in the graphene regions to the Dirac point, V = 
0. The superlattice is contacted at either end by leads at 
potential Vt (not shown). The superlattice potential is plotted 
above the schematic as a function of x. 



In this section, we investigate the shot noise for evanes- 
cent transport through a graphene superlattice consisting 
of many graphene strips tuned to the Dirac point, alter- 
nating with doped regions, as depicted in Fig. [S] We 
eventually take the long superlattice limit of an infinite 
number of such regions. We take the graphene regions 



gated to the Dirac point to be of length L, and they are 
separated by doped regions of length A at the voltage Vd- 
The contacts at both ends are again modeled by highly 
doped graphene at the lead potential Vg. The potential 
Vg drops out of the calculation in the limit Vg,Vd 3> v/L. 
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A. Transmission through a cascade of evanescent 
strips in series 

In order to calculate the transmission probability 
through the graphene superlattice depicted in Fig. [5j we 
employ the transfer matrix method^ The transfer ma- 
trix M(x,x') for the requisite Dirac spinors of a mode 
with transverse momentum q satisfies the equation 



d x M(x,x') = 



V{x) 
-I a x 



M(x,x'). (12) 



In addition, one has^ M(x,x) = I, M(x,x') = 
M(x,x")M(x",x') : detM(x,x') = 1 and 
)a x M(x, x') = a x . The latter condition 
ensures current conservation. 

We write the transfer matrix in the doped regions as 
M d (x,x') = A d {x)A^ l {x l ). Here, A d (x) satisfies Eq. (TT2]) 
in the doped regions, where V(x) — V d . Similarly, the 
transfer matrix in the lead regions contacting either end 
of the superlattice is Me(x,x') = A e (x)Aj 1 (x'), where 
Ag(x) satisfies Eq. (I12p in the lead regions with V(x) — 
Ve. Following Ref. l35l. we choose matrices A e ^(x) that 
allow one to project onto right and left moving states: 



Aj(x) 



where j — {£, d} and k 




-kj+iq c ikiX 



(13) 



with 



Vj/v. 



The transfer matrix in the evanescent regions, 
where the electric potential vanishes, is M e (x,x') — 
exp [qa z (x — x')]. The transfer matrix through a 
single sequence consisting of a region at the Dirac 
point followed by a doped segment is then given by 
Mi(L + A, 0) = M e (L + A, X)M d (X, 0). The transfer ma- 
trix for p doped-evanescent segments in series, contacted 
at either end by ideal leads, is found using Mi and the 
matrices Ag of Eq. (TT31) as 



M p (x, x 1 ) = A e (x)A^[p(L + A)] [Mi(i + A, 0)f x 
x A e (Q)Aj 1 {x') 



(14) 



at x > p(L + A) and x' < 0. 

From the asymptotic form of M p (x,x') in the leads, 
one extracts the transmission amplitudes through the en- 
tire array as explained in Ref. HH: = 1/ \a\ 2 , where 



a P* 
P a* 



lim Aj 1 (x)MJx,-x)AA-x), (15) 

x— >oo 



which, with Eq. (Tl"4"|) . becomes 



a 0* 
P a* 



= A7 x (p(L + A)) [M X {L + A,0)f ^(0). (16) 



The transmission probability T„ ' of mode n is ob- 



tained by setting q 
find 



q n in the above equations, and we 



2A 2 (n) 



[cosh 2 (n£) — l] cosh [2p ■ arcsinh(A(n))] + [cos (2kd(n)X) cosh 2 (n£) — l] 



(17) 



Here, A(n) = y cos 2 [kd(ri)\] cosh 2 (n£) — 1. As previ- 
ously mentioned, the transmission probability becomes 
independent of the lead potential in the limit K^,K d S> 
1/L that we have taken. In the special case of p = 2, this 
equation reduces to Eq. (0. 



The Fano factor can be calculated numerically for a 
cascade of arbitrary length by substitution of the trans- 
mission probability dl~7|) into the integrals 1\ and I2 of 
Eq. (J4j> . The results are plotted as a function of the cas- 
cade size p in Fig. [6] with error bars caused by errors 
in the numerical integration. The Fano factor decreases 
rapidly from its maximum value of 1/3 for p — 1 and ap- 
proaches 0.180... as p — > 00. We next compute the exact 
value of F in this limit p — > 00 of a long superlattice, 
again assuming A >• \n d \ L 2 . 



B. The long superlattice limit 

The transmission probability of Eq. (Ti~7| is plotted 
as a function of mode index n for various values of A and 
p in Fig. [7] One observes a series of peaks in Tn P \ whose 
number increases with A/ | illustrated in the 

first row of plots in Fig. [7] The peaks occur at values of n 
where A is imaginary, and is exponentially damped 
whenever A is a real number. 

When A is imaginary, the second hyperbolic cosine 
in the denominator of Tn^ 1 oscillates with a frequency 
proportional to p, the number of graphene segments in 
the cascade. Therefore, has a series of "subpeaks" 
in the regions of imaginary A, whose number increases 
with p, as illustrated in the second row of plots in Fig. [7] 
We further observe in the second row of Fig. [Jj that the 
damping to the sides of the regions with imaginary A is 
enhanced as p increases. This is also due to the factor of 
p in the second hyperbolic cosine function of Eq. (|17[) . 

Motivated by the above observations, we partition the 
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FIG. 6. Fano factor for a graphene superlattice, as pictured 
in Fig. [5] as a function of p, the length of the cascade. The 
curve is calculated from the transmission probability given in 
Eq. (fTTl) with A = 10\K, d \L 2 . The Fano factor starts at its 
maximum value of 1/3 for a single graphene strip, as con- 
sidered by Tworzydlo et al^ In the long superlattice limit, 
p — > oo, the Fano factor approaches 0.180..., in agreement 
with our analytic calculations. The error bars are caused by 
errors in the numerical integration. 




FIG. 8. Transmission probability Eq. (|17l) as a function of 
mode index n for p — 4 and A = 4 \nd\ L 2 , showing peak os- 
cillations with period set by A and sub-peak oscillations with 
period set by pX. The first peak region to the right of the cen- 
tral region, indexed by n\, is shown magnified in the inset. 
The peak period is given by An\(ri\) in the neighborhood of 
peak n\, while the period of sub-peak oscillation is An p (n p ), 
around peaks indexed by n p . Both periods become infinitesi- 
mal in the limit p ^> 1 with A S> \nd\ L 2 . 
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FIG. 7. Plots of the transmission probability Eq. (|17[) as 
a function of mode index n. The cascade size is p = 2 in 
the first row of plots, while the doped region thickness takes 
the increasing values A = 0.01 |k<j! L 2 (a), A = 2\Kd\L 2 (b) 
and A = 20 \nd\ L 2 (c). Larger thicknesses lead to a higher 
frequency of peaks in the transmission probability. In the 
second row, the doped region thickness is A = 2 \na \ L 2 , and 
the increasing cascade sizes p = 2 (a), p — 6 (b), and p — 10 
(c) cause an increasing frequency of the sub-peak oscillations 
and strong damping outside of the peak regions. 



wave numbers into a series of peak regions where A is 
imaginary. The length of these regions is on the or- 
der of An(n\) in Eq. ([7]). Each peak region is then 
sub-partitioned into sub-peak regions, with a length that 
corresponds to the oscillation period of the second hyper- 
bolic cosine of Eq. (fTTj) at imaginary A. That length is of 
the order An p (n p ) ~ An(n p )/p. Fig. [S] illustrates both 
periods of oscillation, showing the transmission proba- 
bility Eq. (fT7| for p = 4 with doped region thickness 
A-4ML 2 . 

In order to calculate the Fano factor for a long 
graphene superlattice (p ^> 1), we first integrate I\ and 
1-2 of Eq. (|4]) for each peak region An\{n\) over all sub- 
peak regions An p (n p ). The result is then summed over 
the peak region An\(n\) by a second integration, result- 
ing in the desired integral over one peak region. Finally, 
in the limit A 3> L 2 , the peak spacing becomes small 
compared to the decay scale of the transmission proba- 
bility, set by C, and the full integral can be found by fur- 
ther summing the single peak result over all such peaks, 
amounting to a third integration. 

All of the above integration steps can be performed 
analytically, giving the exact Fano factor in the limit 
p -> oo with X/\n d \L 2 -> oo: F = 1/(8 In 2). This re- 
sult is in agreement with the numerical results plotted in 
Fig. [(J and it differs from the Fano factor calculated for 
the geometry of section [Til as wei l as the calculations of 
Tworzydlo et al£ The shot noise for transport through 
a graphene superlattice as depicted in Fig. [S] thus pro- 
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vides another unambiguous signature of evanescent mode 
transport. 



IV. DISCUSSION AND CONCLUSIONS 

Throughout this article, we have assumed translational 
invariance in the y-direction. We now briefly discuss the 
sensitivity of our results to a breaking of that symmetry 
before summarizing. As a typical mechanism for such 
symmetry breaking, we consider a gate edge which makes 
a non-zero angle A0 with the vertical. At such an inter- 
face, transverse modes are mixed. Electrons with mo- 
mentum k x ~ Hd — Vd/v in the ^-direction increase their 
momentum in the y-direction by an amount Ag w Kd^4> 
as they traverse such an edge. 

To avoid qualitatively changing the above calculations, 
this momentum shift must be small on the scale Kd/pXq 
of the oscillations of the transmission probability, requir- 
ing A(j) <C 1 /p\q. Noting that the typical momenta con- 
tributing to the Fano factor in our calculation are of order 
q ~ 1/L, this results in the condition 

A<p < L/pX. (18) 

In addition, we have assumed A 3> \nd\L 2 in all of our 



analytical calculations. These results therefore require 



1<<Mi<< I << ^ 



(19) 



in setups with angular imperfections. 

In conclusion, we have proposed several experimental 
geometries for which the shot noise provides an unam- 
biguous signature of transport through evanescent modes 
at the Dirac point of graphene. For the case of two 
gated graphene strips in series, separated by a long highly 
doped region, the Fano factor can be controlled through 
either spatial asymmetry or the gate potentials. In par- 
ticular, the Fano factor has a universal minimum of 
F = 1/4 for the spatially symmetric case with both gated 
regions tuned to the Dirac point. This result differs from 
the value F = 1/3 expected for transport through a dis- 
ordered metal, which allows for a particularly conclusive 
experimental investigation of evanescent mode transport 
in graphene. For the case of a long superlattice of evanes- 
cent regions we predict another universal value of the 
Fano factor, F = 1/(8 In 2). 
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